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Abstract : Using finite-size-scaling methods, we study the quantum chain version of the 
spin-l-Blume-Capel model coupled to an imaginary field. The aim is to realize higher order non- 
unitary conformal field theories in a simple Ising-type spin model. We find that the first ground- 
state level crossing in the high-temperature phase leads to a second-order phase transition of the 
Yang-Lee universality class (central charge c = —22/5). The Yang-Lee transition region ends at 
a line of a new type of tricriticality, where the three lowest energy levels become degenerate. The 
analysis of the spectrum at two points on this line gives good evidence that this line belongs to 
the universality class of the .A^j-conformal theory with c = —68/7. 
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1 Introduction 



The study of statistical models coupled to imaginary fields can give useful insight on the singu- 
larity structure of the partition functions and the corresponding phase transition mechanisms. 
In 1952 Yang and Lee [|J initiated these investigations considering the Ising model in a complex 
magnetic field B. In the high-temperature regime above the critical temperature T c , the authors 
of fij found that the partition function has zeros for purely imaginary fields B = ih only. In 
the thermodynamical limit these zeros become dense and produce a singularity at h = ±h c (T), 
which for T — > T c gives rise to the Ising phase transition. Fisher Q pointed out that the singu- 
larity in the £>-plane discovered by Yang and Lee has many properties of a standard second order 
phase transition. With the advent of the conformal field theory (CFT)-interpretation of two- 
dimensional phase transitions ||, Cardy Q showed that the Yang-Lee singularity is described 
by the CFT with the simplest possible operator content (just one primary field with conformal 
dimension d = — | apart from the unit operator with d = 0) which has the Virasoro central 
charge c = — -j-. This CFT is non- unitary and corresponds to the case p = 2, p' = 5 of the 
minimal series M PtP i, with central charge 

c = l-6(p-p') 2 /(pp'). (1) 

Generally, a minimal theory M p _ p i (p coprime p') is non-unitary if \p — p'\ > 1, since according 
to the Kac- formula for the scaling dimensions of the primary fields h r s 



h r ,s = V-np-. = ^ 1 - ■ / I' - 1 : 1 - s - ' • s - tx {2) 



it contains at least one primary field with negative dimension h rs . If one wants to describe non- 
unitary minimal models by a coset construction, Kac-Moody algebras at fractional level have to 
be used [||, ^J. For non-unitary minimal theories the modular invariant partition functions and 
their relation to the fusion rules have been studied in M, |8], [| . Recently, such non-unitary models 



appeared, coupled to 2-dimensional gravity, in the study of multicritical matrix models [10|. For 
applications of non-unitary evolution operators in other fields of physics, see e.g. [ |11| . 

For p = 2 the coprime values for p' are p' = 2n + 3, n = 1, 2, . . . (the first member of this series 
is the Yang-Lee-singularity) , while for p = 3 there are two series p' = 3n + 1 and p' = 3n + 2, etc. 
So, the simplest non-unitary minimal theories beyond the -M2,5 are the M.2j which has c = — y 
and A^3 : 5 with c = — |. Lassig ]l2| has considered the renormalization flow between non- unitary 
theories caused by the ^>( 13 )-perturbation. He found that the flow leading to the .M2,5-theory 
is . . . — > ^8,11 -^5,8 - ► -M-2,5- So in connection with the flow to the Al2,5-theory it may be 
useful to check the possibility of a Ai^^ showing up. 

The aim of this paper is to find simple SU (2)-spin quantum chain realisations of non-unitary 
minimal theories. Except for M.2,5, which is well-known to be represented by the Yang-Lee 
singularity of the Ising chain in an imaginary field, no such quantum chain hamiltonians real- 
izing non-unitary A4 pp i theories have been written down in the literature. There is a general, 
though not simple, way for obtaining a quantum chain model for a .Mpy-theory which is even 
integrable: The Forrester-Baxter |l3| 2-dimensional RSOS-models have critical points realizing 



TWpy-theories [14]. Using e.g. the recipes of DeVega [15] one may deduce a quantum chain from 
the RSOS-transfer matrix. However, this has not been done and certainly will lead to very com- 
plicated hamiltonians. In order to get a simple realization of a higher-order non-unitary phase 
transition we shall instead try to generalize the non-integrable Ising- Yang-Lee model, looking 
into its spin-1 analog, and there to possible tripel-level crossings, as will be explained below. 

Let us briefly review the main features of the Ising representation of the .M2,5- m odel. Already 
the work of [Q] and M suggested that the c = — ^ CFT should be approximated on the lattice 



by the thermodynamic limit of an Ising-type quantum chain with the hamiltonian 



1 N 

n LY = -■= J2( a i + X °t° z i+i + ihaf). (3) 
1 i=i 

at a certain curve X(h). Here af and af are standard Pauli matrices acting at site i. X ~ T _1 
is the inverse temperature and h is the imaginary part of the external field. For h = the 
model (H) has the ^-symmetry af <-> —of, af <-> erf, which for h ^ survives as an antiunitary 
symmetry |l6|, 17 1. Due to this, the spectrum of (|3|) remains unchanged under the reflection 
h <-> —h. The operator content and central charge at a particular critical point of (||) were first 
checked by Itzykson et al. [JjJ and were found to agree with the spectrum of the c = — ^ CFT. 
In Hi the model © was studied with high precision by finite-size scaling methods (FFS) (see 



[18] for references to earlier work which used the renormalization group and other techniques). 
Fig. 1 shows the phase diagram of the model (|3|) as it was obtained in [|18|]. 

Although (||) is non-hermitian, it has the special property that the coefficients of its character- 
istic polynomial are real. This is because the antiunitary property of TLly makes it a symmetrical 
matrix. Therefore the eigenvalues of © are either real or come as complex conjugate pairs. 

For h = 0, TLly is hermitian and the spectrum is real. In the high-temperature regime 
A < 1 there is a finite gap. By continuity, in this regime, the finite- N energy levels can become 
complex only if they cross and pairwise form complex conjugate pairs. This happens first just at 



the c = — phase transition and this property was used in [18| to obtain the FFS determination 



5 

of the phase transition curve shown in Fig. 1. For h = and A = 1 we have the standard Ising 
phase transition described by the minimal CFT with c = I. For h = and A < 1 we are in the 
massive ^-symmetric high-temperature phase. 

Within a mean-field theory or Landau approach the appearance of a phase transition for 
A < 1 and h ^ can be explained as follows 0: Consider a model for which the free energy 
admits a Landau expansion in terms of an order parameter tp which couples to a magnetic field 
B: 

F = Fo+Bv + V 2 f2 + + • • • (4) 

J^o, T% and JF4 are functions of the inverse temperature A. For B = a second order phase 
transition will take place if there is a A = A c with J-2(X C ) = and J-i(X c ) > 0. Now, for any 
A = Xly with ^(Xly) > and ^(Xly) > (this is in the high-temperature phase) criticality 
of the system can be restored by applying a correctly tuned imaginary magnetic field: shifting 
the order parameter by the constant iipo, i.e. (p (p + ipo, (0) becomes 

T = B{(p + icpo) + (<p + i(po) 2 T2 + (<p + ^o) 4 ^4 + • • • 

= ... + p[B + 2ip Q {T 2 - 2p Q T A )] + $ 2 {T 2 - 6^4) + 4^Vo^4 + -.., (5) 

where the dots stand for terms of order ip°, and higher than (p 3 . The coefficient of ip 2 will vanish 
if we choose 



P0 = ± S JJ 7 2{X L Y)/m{XLY), (6) 

giving rise to the Yang-Lee phase transition. This way, necessarily a term proportional i^ 3 is 
introduced. Neglecting the residiual </? 4 -term, at A = A^y we have 



T= (B-B L y)p + ±igp 3 + ... (7) 



with 



3 



BLY = ~Jym and 9 = ±V24^4. (8) 



We now address the question, whether and how can we get non-unitary multicritical phase 
transitions in the sense that not only the coefficient of (p 2 but also higher terms of a Landau 
expansion come to vanish. Since, in order to achieve this, we will need more degrees of freedom 
than are available in a spin-^ chain, we shall study a spin-1 quantum chain placed in an imaginary 
field which breaks the original ^-symmetry. Is there a new type of phase transitions which may 
be called a non-unitary tricritical phase transition and which is produced by the crossing of the 
three lowest energy levels? In general we should have three parameters in the hamiltonian in 
order to find such a triple cross-over. Looking into the Fisher's Landau expansion @) this has 
a chance to be realized in a model in which the parameters can be choosen to get vanishing 
coefficients of if 2 , if 4 and also of if 3 , so that the first nonvanishing terms is proportional c^ 5 (with 
a purely imaginary coefficient). In a mean field treatment of the non-hermitian spin-1 model 
defined in (^) below, this tricritical curve has been determined by K. and M.Becker [[!(], [lTj. So 
it seems possible to generalize the realizations of the (A m -i,A m ) unitary minimal CFTs (which 
have \p — p'\ = 1, p' = <m > 3), by spin-s-quantum chains with s = m/2 — 1 [16, 17, |l9| to the 
non-unitary series. 

Attempts to establish a connection between the mean-field non-hermitian multicritical curve 
and a particular minimal CFT via a generalization of Zamolodchikov's composite field repre- 
sentation of the CFT fusion rules, have met various difficulties [[l(], 17]. So, as a more direct 
approach, in this paper we treat a suitable non-hermitian quantum chain hamiltonian by FSS 
methods and determine its non-hermitian tricritical manifold. From the chain size-dependence 
of the low-lying levels of the spectrum we then determine the CFT's which describe the spectrum 
in the different regions of the critical manifold. 

The main part of this paper is organized as follows: In Sec. 2 we define the particular spin-1 
quantum chain we want to study and review its known zero magnetic field properties. In Sec. 3 
we determine the phase transition in the imaginary field from the ground-state level crossings. 
Sec. 4 is devoted to the calculation of the finite-size spectrum on the two-level crossing phase 
transition curve and the determination of the universality class involved. Sec. 5 contains the main 
result of the paper: the determination of the FSS of the spectrum at the tripel level crossings 
which gives the evidence that there the .M2,7-minimal theory is realized. Sec. 6 presents our 
Conclusions. 



2 The spin-1 Blume-Capel quantum chain with an imaginary 
field 

We consider a spin-1 quantum chain which has the property that the characteristic polynomial 
of the hamiltonian matrix has only real coefficients, so that also in this case the eigenvalues are 
either real or occur in complex conjugate pairs. This property seems necessary if we want to 
find universality classes described by minimal CFTs, because there the scaling dimensions of the 
fields, which are related to the energy gaps of the quantum chain hamiltonian, are always real. 
We insist in having a ^-symmetry for non-zero imaginary magnetic field, because the simplest 
modular invariant partition functions, those of the p = 2-series, exhibit a ^-symmetry. 
The non-hermitian symmetric spin-1 hamiltonian we shall consider is defined by 

N 

H = Y, («(^) 2 + PSf - S?Sf +1 - ihSf) . (9) 
i=i 

Here Sf and Sf are standard 3x3 spin-1 matrices acting at site i. We shall choose Sf to be 
diagonal. This 7i has an antiunitary ^-symmetry: Denoting by K the operator of complex 



conjugation, and denning [16 







A? 



u = n (ks x ? - 1) , (io) 

i=l 1 

TC satisfies 

[n, e] = o with e = ku. (11) 

There are three parameters: a, /3 and /i. For h = this is the standard Blume-Capel [^] 
quantum chain, and certain linear combinations of a and (3 can be interpreted as a temperature 
variable and a vacancy density. In principle, one may add more ^-symmetrical terms to (||), 
e.g. terms proportional to (Sf) 2 and to (Sf ) 2 (Sf +1 ) 2 etc. (Blume-Emery-Griffiths model pl[). 
However, for simplicity, here we shall consider the form (|9|) only. 

Let us first review the phase diagram of (|9|) for h = 0. The critical behaviour of the Blume- 
Capel quantum chain (with zero field h) has been studied first in mean field approximation by 
Gefen et al. pO] and then, using finite-size scaling (FSS) by Alcaraz et al. p^, |23| and by one 



of us 24], see also [25]. The line (3 = is classical and is easily seen to contain a first order 
transition at a = 1. The spectrum of (0) is unchanged by the reflection (3 — > —j3. So it is 
sufficient to consider only (3 > 0. Using the FSS-technique, the transition line for [3 ^ can be 
located quite precisely as can be seen from Table 1. 



a 


-1.0 


-0.5 


0.0 


0.20 


0.50 


0.55 


0.60 


0.65 


0.70 


0.725 


Pc 


1.884 


1.6266 


1.3259 


1.187 


0.9411 


0.895 


0.8448 


0.7932 


0.737 


0.707 


a 


0.76 


0.80 


0.84 


0.85 


0.90 


0.9102 


0.93 


0.95 


0.98 


1.00 




0.6618 


0.6070 


0.5458 


0.531 


0.4365 


0.4157 


0.368 


0.3118 


0.199 


0.0 



Table 1: Critical values of the Blume-Capel-quantum chain, eq.(||) with h = 0, as obtained 
from finite-size scaling |24[| . The underlined values of a,/3 are those of the tricritical point. 



There is a single critical curve which starts at a = 1 and (3 = and then moves with increasing 
(3 towards smaller values of a. It first follows the unit circle in the /3-a-plane and then takes off 
to large negative a, ultimately following the curve (3 = \/2\a\. Still close to a = 1, it passes a 
tricritical point at [24]] 



a t = 0.910207(4), (3 t = 0.415685(6). (12) 

For a < at the phase transition is second order, and, on the transition line, the low-lying part 
of the spectrum is that of the c = \ Ising CFT. The high-lying spectrum should also contain 
a massive component, since the Ising Maiorana fermion needs only two of the three degrees of 
freedom which are available. At the tricritical point the spectrum is described by the c = jq 
CFT ("tricritical Ising model"). For at < a < 1 the transition is first order. The phase below 
the critical curve at a < 1 and (3 > is a low-temperature phase with broken .^-symmetry. 



3 Determination of the phase transition curves from lowest level 
crossings 

In analogy to the situation for the Ising chain @, we expect that also in the high-temperature 
regime of (^) there should be a range in h where the spectrum remains real despite the hamiltonian 

^The error in the last given digit of a number will be quoted in brackets, e.g. 0.910207(4) stands for 0.910207 ± 
0.000004, etc. 



being non-hermitian. The spectrum should become complex where the first lowest level crossings 
appear. 

We have calculated the four lowest energy levels of the translationally invariant states of the 
hamiltonian @ with N = 2, . . . , 11 sites, for various values of the parameters a, (3 and h, using a 
Lanczos technique. For JV<8we diagonalize the hamiltonian also by exact methods without the 
restriction to momentum-zero states. In this paper we always take periodic boundary conditions. 

As a typical example of the lowest level crossing behaviour at our small values of N, in the 
last line of Table 2 we show results for a fixed value of a, i.e. a = 0.80, and the value (3 = 0.725, 
which for h = would be in the disordered phase, slightly above the critical curve (for h = 
from Table 1 we have (3 C = 0.6070). 

Two-level cross-over 

Now, switching on h, we find that as long as h ~ 0.24, the four lowest energy levels remain real as 
they were for h = 0. However, e.g. for N = 10 sites at h c = 0.025828 the two lowest levels meet 
and form a square-root branch point. For h < h c the two eigenvalues of the spectrum with the 
lowest real parts form a complex conjugate pair. The branch point positions are iV-dependent. 



p 


a = 0.80 

7V = 5 6 7 8 9 10 11 


oo 


0.675 
0.700 
0.725 


0.023796 0.019666 0.016984 0.015119 0.013758 0.012725 0.011920 
0.029975 0.025677 0.022893 0.020961 0.019553 0.018486 0.017655 
0.037034 0.032733 0.030000 0.028142 0.026815 0.025828 0.025072 


0.007(1) 
0.012(1) 
0.021(1) 





a = 0.20 

iV = 56 7 8 9 10 11 


oo 


1.800 
6.000 


0.143269 0.132717 0.126115 0.121739 0.118711 0.116539 0.114938 
2.141181 2.124866 2.115361 2.109411 2.105475 2.102753 2.100700 


0.1080(2) 
2.0931(2) 



Table 2 : Values h c (N) of the first crossing of the two lowest levels of the hamiltonian @ for 
a = 0.80 and three values of (3, and for a = 0.20 and two values of (3. For h above these values 
the two lowest eigenvalues form a complex conjugate pair. N denotes the chain length, the last 
column N = oo gives an estimate of branch point position in the thermodynamic limit. 



We use standard extrapolation techniques (both rational polynomials and Van-den-Broeck- 
Schwarz approximants f26f ) to obtain the branch point positions in the limit N — > oo. In 
analogy to the situation in the Yang-Lee- model @ [18], we interpret these as second-order phase 
transition points. 

Table 3 collects our numerical results for these branch point positions h c (a,(3) obtained by 
extrapolation h c = lmi7v->oo h c (N) from N = 3, ... ,11 sites. In order to save space, in Table 2 
we give some abbreviated data, but for the actual extrapolation we used also iV = 3 and N = 4 
sites and 8-digit precision for the finite- iV-branch point positions. 

In Fig. 2 we plot the resulting critical surface in the space of a, (3 and h. This surface has 
the shape of wings starting from the h = 0- critical curve of Table 1. The wings converge at a 
vertex which is the standard tricritical Ising point (p^) with c = jq. This shape of the wings 
is just the same as it is familiar for applied staggered real fields, see e.g. Fig. 44 in Lawrie and 
Sarbach lEl. 



a = 0.20 

P h c 


a = 0.50 

P h c 


a = 0.55 

P K 


a = 0.65 

P K 


a = 0.725 

P h c 


a = 0.76 

P K 


1.187 0.0 
1.8 0.108(1) 
3.0 0.561(1) 
6.0 2.08(1) 
9.0 3.81(1) 
18. 9.40(2) 


0.941 0.0 
1.8 0.2172(2) 
3.0 0.7558(2) 
6.0 2.409(2) 
9.0 4.21(1) 
18. 9.96(1) 


0.895 0.0 
0.95 0.0025(6) 

1.0 0.0070(5) 

1.1 0.0228(6) 

1.2 0.0445(4) 
1.8 0.246(3) 
3.0 0.807(3) 


0.792 0.0 
0.85 0.0027(3) 
0.95 0.0169(8) 
1.0 0.0273(6) 
1.2 0.0820(4) 
1.5 0.1922(6) 
1.9 0.374(2) 


0.706 0.0 
0.75 0.0020(2) 
0.80 0.0085(3) 
0.85 0.0176(5) 
0.90 0.030(1) 
1.0 0.059(1) 


0.662 0.0 
0.70 0.0018(2) 
0.75 0.0088(3) 
0.80 0.0194(3) 
0.85 0.033(2) 
0.90 0.050(1) 



Table 3 : Non-hermitian critical values h c (a,P) obtained from the lowest ground state level 
crossings, except for the first line, which gives the Blume-Capel value at h c = 0. 

Triple cross-over: Tricritical line 

In calculating the ground-state level crossings of the Yang-Lee model @ we don't have to worry 
about the third and higher energy levels, because there always the first two levels meet before a 
third level crosses in. In the spin-l-case the situation is different because we are varying one more 
parameter. For fixed 0.55 < a < at, moving out in p, we reach a point, which we shall call Pt, 
where the third energy level E% crosses the second level E\ at a value of h which is lower than the 
crossing point of the two lowest levels Eq and E\. In Fig. 3 we show this situation for N = 6 and 
a = 0.8, where it occurs for P ~ 0.78097. For a fixed finite chain size N, there is no point where 
all three lowest levels meet: the tripel crossing is an avoided cross-over since in our case there is 
no new conservation law at this special point. Table 4 gives our determination of the positions 
of these avoided cross-overs at chains of N sites and several values of a. The calculation for 
obtaining these data requires a considerable amount of computing time and therefore we mostly 
have not pushed the precision beyond four significant figures, although this would be possible 
if wanted. The value of ht quoted is taken somewhat arbitrarily as the lowest value of h where 
at P = Pt we have E\ = E%, see Fig. 3. The convergence of the values of Pt and ht with N is 
quite slow. So there is a considerable uncertainty for the values extrapolated to N — > oo. It may 
be that for a ~ 0.5 both Pt and ht move to infinity, so that for low a the wings may extend to 
infinity. 

For a — > at the values of ht go to zero. Accordingly the width of the wings (for N — ► oo) is 
seen to shrink towards zero. This means that the range in h where the spectrum remains real 
decreases if a moves towards the tricritical point. This is natural, since in the adjacent lst-order 
region there is a ground state degeneracy in the limit N — > oo at h = 0. 



a 


N 


3 


4 


5 


6 


7 


8 


9 


10 


00 


0.55 


Ht 


8.8958 


7.63 


7.45 


7.53 


7.66 


7.82 


7.887 


7.98 


? 




h 


44793 


3.61 


3.46 


3.49 


3.56 


3.622 


3.681 


3.73 




0.60 




4.2604 


3.64 


3.462 


3.409 


3.400 


3.4022 


3.413 


3.40 


? 






1.7276 


1.312 


1.185 


1.140 


1.123 


1.198 


1.121 


1.11 




0.65 




2.602 


2.2336 


2.1035 


2.0485 


2.0235 


2.0123 


2.0075 


2.0060 


2.006(6) 




h t 


0.852 


0.613 


0.5249 


0.484 


0.464 


0.4535 


0.4474 


O.444 


0.437(4) 


0.725 


Pt 


1.499 


1.3094 


1.2310 


1.191 


1.1693 


1.156 


1.1480 


1.1428 


1.127(3) 




h t 


0.356 


0.2363 


0.187 


0.161 


0.146 


0.137 


0.1309 


0.1268 


0.113(2) 


0.76 


0t 


1.204 


1.062 


1.005 


0.9697 


0.951 


0.9391 


0.9315 


0.92596 


0.908(5) 




h t 


0.246 


0.1566 


0.119 


0.0992 


0.0872 


0.0796 


0.0743 


0.07062 


0.060(3) 


0.80 


A 


0.9482 


0.849 


0.806 


0.780 


0.7663 


0.7566 


0.7501 


0.74524 


0.728(7) 




h t 


0.1601 


0.098 


0.071 


0.0569 


0.048 


0.0419 


0.0382 


0.03529 


0.025(3) 


0.84 


Pt 


0.748 


0.6814 


0.651 


0.636 


0.624 


0.617 


0.6120 


0.6083 


? 




h t 


0.103 


0.060 


0.0415 


0.031 


0.025 


0.0217 


0.0188 


0.0167 





Table 4: Finite- ^-approximations to the non-hermitian " tricritical" points of the model (||). 
In brackets we indicate the estimated error in the last given digit. If no bracket is given, then 
the error is of the order of one unit in the last given digit. In the last column "00" we give the 
extrapolated value for N — > 00 for those cases in which we get reasonable convergence. 



In Fig. 4 we show the projection of the critical curves for fixed values of a in the — /i-plane. 
The curve connecting the end points is the tricritical curve. In the mean field treatment of (||), 
which was mentioned earlier in the Introduction, in the neighbourhood of the tricritical point 
K. and M.Becker [16, [l7| have obtained curves very similar to those of Fig. 4. They obtain the 
endpoints of the Yang-Lee-transition curves through the vanishing of the Landau-coefficient of 
ip 3 in the notation of (|5|), (R). Since in mean- field the tricritical point for the h = model is at 
(a, 0) = (0.5836 . . . , 1.2011 . . .) instead of the correct value from FSS (0.91021 . . . , 0.41569 . . .), 
their corresponding curves are at too small values of a and about three times too large values of 
0. Anyway, it is interesting that near to the tricritical point, mean field reproduces the shape of 
the fixed-a critical curves so well. 



4 Determination of the spectrum on the phase transition surface 

In order to determine the universality class of the non-hermitian phase transitions found, we use 
standard finite-size-scaling (FSS) methods of conformal theory. Using always periodic boundary 
conditions, from the iV-dependence of the ground-state energy E${N): 

E (N)/H = -Na - — + ... (13) 

we obtain the effective central charge c, which for non-unitary theories is given by 

c = c-24h min , (14) 

with c the central charge and h m % n = (1 — (p— p') 2 ) / {^pp') is the lowest (most negative) anomalous 
dimension, ao is the bulk constant which will not be of interest in the following, £ is the conformal 
normalization factor 28]. Inserting h m i n into (|l|) we can write 



2=1 — 6/(pp ). 



(15) 



We calculate N sites approximants Cjy to c from 



CAT 



6N(N — 1) 
(2N - 1)tt£ 



N - I)Eq(N) - NE {N - 1 



(16) 



and extrapolate N — ► oo. At a conformal point the gaps scale in leading order of N as N . The 
coefficient of N~ 1 allows us to obtain the anomalous dimensions (h, h) of the primary conformal 



fields |28fl: 



Xi{P) 



lim Xi(N, P) = lim g t (N, P)/£ = {h + r + h + r)< - (h + h) r 



with 



N 

Gi(N, P) = —(Ei(N, P) - E (N)). 

Z7T 



(17) 



(18) 



Here Ei (N, P) denotes the energy of the i-th level in the momentum P-sector of the hamiltonian 
with N sites. The positive integers r, r become nonzero for descendant field levels and determine 
the momentum through P = r — r' . The term (h + h) m i n appears in case of non-unitary theories, 
since there the ground state does not correspond to the vacuum. 

The conformal normalization £ at the critical point under consideration is calculated from 
either the lowest AP = 2-gap in the thermal sector, or from the lowest AP = 1-gap in the 
non-zero-charge sector [28]. If possible, we usually use both gaps in order to have a consistency 
check. 

Tables 5 and 6 give our finite- iV-results for the scaled gaps of the lowest levels at two typical 
values of the critical surface below ht (i.e. on the wings) (a = 0.2, (3 = 6.0, h = 2.0931) and 
(a = 0.5, (3 = 3.0, h = 0.7559), together with the values expected for c = — We see that we 
have excellent agreement with the c = — ^ conformal spectrum. Our spin-1 model has one more 
degree of freedom than the spin-^-Ising chain, but probably the extra degree of freedom cannot 
be seen in our data because the corresponding levels are massive and are high above the ground 
state. 
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Table 5: Finite-size data for the scaled gaps x(N,P) for a = 0.2, (5 = 6.0, h = 2.0931 and 
the normalization £ = 2.84. In the bottom line we give the values expected for the universality 
class is M.2 5- 









p — n 

/ — u 






P — 1 




P — 9 


/V 

i V 


C 


X\ 


x 2 


~ 

x 3 




^2 


x 3 


X\ 


9 




U.ZUOc/Uc/ 


U.UO / li 


1 9944 


u.uuuuu 


U.OOU 1 






Q 
O 


D ^1 79RQ 


u.ou i yyu 


D 89^88 


1 ^988 

l.tJiOO 


fl 8fi7zL8 


1.UU1U 


1 R88R 


u.ou / o 


4 


f) 463568 


336676 


98033 


1.7606 


94030 


1.2531 


2.1169 


1.3479 


5 


0.440614 


355312 


1.12429 


1.8785 


97035 


1.4288 


2.3901 


1.5953 


6 


0.429010 


0.367693 


1.25671 


1.9671 


0.98507 


1.5934 


2.6054 


1.7750 


7 


0.422552 


0.376158 


1.37522 


2.0661 


0.99321 


1.7491 


2.6931 


1.9181 


8 


0.418705 


0.382111 


1.47783 


2.1709 


0.99812 


1.8964 


2.7785 


2.1921 


9 


0.416309 


0.386402 


1.56387 


2.2825 










10 


0.414784 


0.389560 


1.63435 


2.4001 










11 


0.413823 


0.391920 


1.69136 


2.5220 










12 


0.413249 


0.393703 


1.73742 


2.6457 










oo 


0.38(1) 


0.42(1) 


2.02(3) 


4.3(1) 


1.011(2) 


3.35(3) 


3.09(6) 


2.6(5) 




0.4 


0.4 


2.0 


4.0 


1.0 


3.0 


2.0 


2.4 



Table 6: Finite-size data for the scaled gaps x(N, P) as in Table 5, but for a = 0.5, (5 = 
3.0, h = 0.7559 and the normalization £ = 1.90. For P = 2 the lowest levels are real for 
N = 3 and N = 8 sites, for N = 4, . . . , 7 they form complex conjugate pairs. This makes the 
extrapolation rather uncertain. 



As another check, whether the wings are a critical surface of the c = —^-universality class, we 
have also looked into the behaviour of the mass gap in the neighbourhood of the critical surface. 
Since the only relevant perturbation of the c = — ^-conformal theory is by the primary field 
_i) (its descendent is redundant) [|], 29], we expect the mass gap to increase proportional 

to r 5 / 12 where r is the distance from the critical curve. The particular choice of the direction in 
which we take this distance is not relevant, since all directions are equivalent in leading order, 
see the analogous situation in the Yang-Lee Ising model [18|. Since the perturbed c = — 
theory contains only one single particle of mass m(a, (3,t), it follows that if the first massive 
level is at AEi = m, the second massive level should appear at AE2 = 2m, followed by a bunch 
of levels close above 2m. Table 7 gives off-critical masses (extrapolated to JV -> 00) for three 
approaches to the critical surface. These data confirm the expected ^-power law and the relation 
AE2 ~ 2AE\ very well. Observe that also in the last example, a = 0.65, (5 = 2.007, h — > 0.431, 
which is at an end point of the critical c = —^-surface, the data clearly show the ^ -power 
behaviour. The error bars are much larger for AE2 since for two-particle states the convergence 
is no longer exponential in N, but only powerlike N~ 2 . Finite lattice effects also limit the validity 
of the ^-power law to the neighbourhood of the critical point [30]. 
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Table 7: The two lowest gaps AE\ and AE2 in the off-critical regime at various off-critical 
distances converging to the respective critical points, together with, in the bottom row, fit 
formulae for AE 1 s=s \AE 2 . 



5 The spectrum at "nonunitary tricriticality" . 

From the triple crossing points determined for iV = 3 ... 10 sites, see Table 4, in a few cases we 
were able to give an estimate of the limiting values for N — > 00 (see the last column of Table 4) . 
We now look into the spectra at two of these values in order to see whether we can identify a 
particular modular invariant or universality class. Taking the central values given in Table 4, we 
find no good evidence for conformal invariance in the sense that for N = 2, . . . , 12, the NAEi are 
not clearly seen to tend to constants. In order to improve the determination of the triple points, 
we did the calculation of the spectra for a net of about ten points around the central values found 
in the last column of Table 4. We looked at the respective A-dependence of the effective central 
charge cat and of the NAEi and selected the value of a, (3 and h with best convergence in N. 
This way, we were able to improve the determination of the triple crossing points for A" — ► 00 
such that the A^-behaviour required by conformal invariance is seen. 

More precisely, we have calculated the spectra for a = 0.65 at 11 points of the grid: f3 = 
2.003; 2.007; 2.008; 2.010 vs. h = 0.428; 0.431; 0.432; 0.438 and choose the best convergence 
points. Actually, it is not worthwhile to calculate the spectra at all points of the grid, because in 
the corners, convergence is deteriorating fast, so that there we certainly move away from the con- 
formal point. Analogously, for a = 0.725 we calculated at j3 = 0.105; 0.108; 0.110; 0.114; 0.122 
vs. h = 1.126; 1.127; 1.128. Tables 8 and 9 give the finite-size data at the best convergence 
points. 
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Table 8: Finite-size spectra at the extrapolated triple crossing point a = 0.65, (3 = 2.007, 
h = 0.432, using for the normalization factor the estimate £ = 2.1. The imaginary parts are 
given as Ji = ImEi(N,p)/^, i.e. not multiplied by the size factor N/2n. For P = 0, xi,X2 and 
X5 are real, but X3 and X4 come as a complex conjugate pair. For P = 1, the four lowest levels 
come in complex conjugate pairs, while X5 is real. In the line N = 00 we give estimates for the 
thermodynamic limit, which, however, for c and x.2 are so uncertain that we are unable to give 
any reliable estimate. In the bottom line we list the values expected for the universality class 
M 2 ,7- 
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Table 9: Finite-size spectra as in Table 8, here for the extrapolated triple crossing point 
a = 0.725, (3 = 1.126, h = 0.110. We use the estimated normalization factor £ = 1.2. Only for 
N > 7 we have Rex4 = Rex%. As in Table 8, we find that for P = the x±, X2 and X5 are real. 



We notice that there is no general tendency for the Xj or c to converge towards zero, as 
it would appear if we missed the critical point considerably, and instead were in the massive 
regime. However, unfortunately, the sequences don't converge very well in the range of sites 
(N = 2, ... , 12) accessible to our calculation, but one can see a clear pattern of levels emerging. 



We shall now compare these patterns observed in Tables 8 and 9 to the anomalous dimensions 
of A^2,5 and three next simple minimal non-unitary field theories -M p y beyond the M.2,5- 
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In order not to get lost in too many or somewhat exotic possibilities, we shall not consider non- 
minimal models. We shall consider only (A p _i, A p '_i)-modular invariants so that x\ = 2/ij. This 
will be justified by the success of getting along with a minimal theory. From the Kac- formula 
for the anomalous dimensions (§) we calculate the conformal grids: 
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Table 10: Conformal grids of the anomalous dimensions h r 
and A^5,8 (in the grid for A4^ t s we have omitted the row s 
row for r = 1). 



for the theories M.2,5, M.2,7, -^3,5 
7 which is equal to the reversed 



The patterns of the spectra expected for the three cases are quite different, and different 
too from the spectrum of the Yang-Lee theory which we found on the critical wings. For 

the sequences c; x\, £2, £3, ... in the P = and P = 1-sectors we expect (we indicate the 
vacuum levels by underlined numbers): 

. M 2 , 5 : §; f, 2, 4, 4 + §, 6, 6 + §, 8, ... 

c = 0.4; P = : 04, 2.0, 4.0, 4.4, 6.0, 6.4, 8.0, ... 
P =1 : 1.0, 3.0,... 

•M 2J : f; f, f, 2, 2 + f, 4, 4+f, 4+f, 6, 6+f, ... 

c» 0.571; P = : 0.286, 0.857 . 2.0, 2.286, 4.0, 4.286, 4.857, ... 
P = 1 : 1.0, 1.286, 3.0, 3.286, ... 

• ^3,5 : f; jjj, 2> I' 2 ' 2 + i 2 +f' ••• 

c = 0.6; P = : 01, 0.5, 1.6, 2.0, 2.5, 3.6, 4.0, 4.1, ... 
P =1 : 1.0, 1.5, 2.6, 3.0, ... 

m KA • 11 • JL J- J- X 3 1 3 143 9 oil 

^^15,8- 20' 80' 10' 16' 16' 5' ' 2' 80' ' 80 ' '■■ 

c = 0.85; P = : 0.0375, Ol, 0.1875, 0.4375, 0.6, 1.0, 1.5, 1.7875, 2.0, ... 
P = 1 : 1.0, 1.0375, 1.1875, ... 



A first hint about which modular invariant is realized at the tripel crossing point can be obtained 
by comparing c to x\\ Both are equal for the Yang-Lee case A^2,5- For 7^2,7 we expect c = 2x, 
while we should have x <C c for A4^^ and the more for Ai^^- Our finite-size spectra are 
compatible with the second case. The normalization by Rexf^ 1 is quite stable which encourages 
us to look into the absolute value of x\ and X2- Only for M2J the order of magnitude is 
compatible with the data. 

Also the number of levels well below Xi = 2 is quite different for the universality classes 
considered: one P = 0- level in the case of A42,5, two for A42,7, three for .M3.5 and 7 to 8 for 
Msg. Since we observe two levels in this range, M?,^ is not completely out from this simple 
counting since the level at x = 1.6 might have been misidentified, but we find that all data 
together are best compatible with M.2,7- The complex conjugate pairs in P = then must 
come from the first decendents of the two primary fields with h = | and h = | together. The 
appearance of the imaginary parts probably is just a transitory phenomenon at small values of 
N. A further check of the correct assignment is given by the P = 1-levels, which appear at the 
expected positions, although again in complex conjugate pairs. 

We see that our finite-size data strongly hint that the tripel crossing points are described by 
a conformal field theory with c = — A good further check would be calculate the spectrum 
of the model with Z2-twisted boundary conditions |31, ^] and to look whether it corresponds 
to the twisted modular invariant. It is quite obvious to speculate that at analogous four- fold 
crossings in a higher spin version of (|9|) if it contains one more parameter, the A42.9-CFT may 
be realized. 

The off-critical massive theories obtained from the A42,2n+3-theories by their ^^-perturba- 
tions has been worked out in [33, 34]. The S-matrix contains N particles with Koeberle-Swieca- 
mass ratios [35], so that the 0( lj3 )-perturbed A42,7-theory should have two particles with mass 
ratio milmx = sin(27r/5)/sin(-7r/5). Since we have not tried to determine the directions of the 
perturbations in the a, (3, /i-plane, we cannot confirm or disprove this feature in our spin model. 



6 Conclusions 

Using numerical finite-size scaling techniques, we have shown that the spin-1 Blume-Capel quan- 
tum chain in an imaginary field h, eq.(|9|), has a two-dimensional critical surface, arising from 
ground-state level crossings, which belongs to the Yang-Lee-edge A42,5 (c = — ^-universality 
class. On both sides, this surface ends in a new type of tricritical line due to tripel lowest level 
crossings. Fig. 2 shows the phase diagram as determined by finite-size scaling. 

The FSS-spectra at the non-hermitian tricritical lines show the particular pattern expected 
for the A^j-niodular invariant partition function, from which we conclude that these realize the 
c = — y universality class. This is the first time that in a simple 5C/(2)-spin quantum chain 
critical behaviour corresponding to a CFT with c < — ^ is observed. 
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Figure Captions 

Fig. 1: Phase diagram of the Ising quantum chain in an imaginary field h defined by the 
hamiltonian (^), from |18| . 

Fig. 2: Three-dimensional view of the phase diagram of the hamiltonian (||), as obtained by 
finite-size scaling. The five-cornered star at a ~ 0.910, (3 ~ 0.4157 is the tricritical point with 
central charge c = 7/10. It is the origin of a pair of wings with c = —22/5 which extend to 
imaginary fields h from the Ising-like line at h = and a < 0.910. 

Fig. 3: Behaviour of the real part of the three lowest eigenvalues of the spin-l-hamiltonian 
(|j|) near the tripel cross-over region a = 0.80 and j3 « 0.781, for N = 6 sites (see the en- 
try: 0.780; 0.0569 in Table 4). For j3 < 0.78098, Eq and E\ form a complex conjugate pair near 
h rj 0.0566. If P > 0.78098, then E\ joins Ei to form a complex conjugate pair above h m 0.0568. 

Fig. 4: Plot of the wing critical curves obtained from FSS as in Fig. 2, but now for vari- 
ous fixed values of a projected onto the (3 — h plane. 
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